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Abstract 

Using Rule 126 elementary cellular automaton (EGA) we demonstrate 
that a chaotic discrete system — when enriched with memory - hence ex- 
hibits complex dynamics where such space exploits on an ample universe 
of periodic patterns induced from original information of the ahistorical 
system. First we analyse classic EGA Rule 126 to identify basic charac- 
teristics with mean field theory, basins, and de Bruijn diagrams. In order 
to derive this complex dynamics, we use a kind of memory on Rule 126; 
from here interactions between gliders are studied for detecting stationary 
patterns, glider guns and simulating specific simple computable functions 
produced by glider collisions. 

Keywords: elementary cellular automata, memory, Rule 126, gliders, 
glider guns, filters, chaos and complex dynamics. 
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1 Introduction 



In this paper we are making use of the memory tool to get a complex system 
from a chaotic function in discrete dynamical environments. Such technique 
takes the past history of the system for constructing its present and future: the 
memory [HI [3 [31 HI HI [S] . It was previously reported in [TH] how the chaotic EGA 
Rule 30 is decomposed into a complex system applying memory on this function. 
Recent results show that other chaotic functions (Rule 86 and Rule 101) yield 
complex dynamics selecting a kind of memory, including a controller to obtain 
self-organization by structure reactions and simple computations implemented 
by soliton reactions [T7]. 

We focus this work on cellular automata (CA) evolving in one dimension, in 
particular taking the well-known EGA where each function evaluates a central 
cell and its two nearest neighbors (from left and right) and, every cell takes a 
value from a binary alphabet. Such EGA were introduced by Wolfram and have 
been widely studied in several directions [36 . 

Among such EGA, there is a set of functions evolving in chaotic global 
behaviour where a number of cells remain unordered and their transitions have 
a large number of ancestors [38 . In this sense, special attention is given on a 
chaotic one: the EGA Rule 126 introduced by Wolfram in ^4|. Particularly in 
[37]j^it is commented that Rule 126 generates a regular language with average 
growing faster than any polynomial time. This property can be analysed as 
well by de Bruijn diagrams with additional features; in particular another partial 
analysis with de Bruijn and subset diagrams was done by Mcintosh [23; showing 
that string 010 in Rule 126 is the minimal Garden of Eden configuration (having 
no ancestors). 

Based on the idea developed in previous results for obtaining complex dy- 
namics from chaotic functions selecting memory and working systematically, it 
was suspected that a complex dynamics may emerge in Rule 126 given its re- 
lation to regular languages; making use of gliders coded by regular expressions, 
as it was studied in Rule 110 "26] and Rule 54 pT]. 

In this way Rule 126 provides a special case of how a chaotic behaviour 
can be decomposed selecting a kind of memory into a extraordinary activity of 
gliders, glider guns, still-life structures and a huge number of reactions. Such 
features can be compared to Brain Brian's rule behaviour or Gonway's Life but 
in one dimension; actually none traditional EGA could have a glider dynamics 
comparable to the one revealed in this EGA with memory denoted as 4>Ri2emaj 
(following notation described in [IHIIIT]). 

1.1 One-dimensional cellular automata 

One-dimensional GA is represented by an array of cells Xi where i € Z (integer 
set) and each x takes a value from a finite alphabet S. Thus, a sequence of 
cells {xi\ of finite length n describes a string or global configuration c on S. 

^http: / /mathworld. wolfram.com/Rulel26. html 
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This way, the set of finite configurations will be expressed as S". An evolution 
is comprised by a sequence of configurations {ci} produced by the mapping 
$ : — >■ S"; thus the global relation is symbolized as: 

$(c*) ^ c*+i (1) 

where t represents time and every global state of c is defined by a sequence of 
cell states. The global relation is determined over the cell states in configuration 
c* updated at the next configuration c*^^ simultaneously by a local function ip 
as follows: 

(p{xl_,,...,xl...,xl+,)^xl+\ (2) 

Wolfram represents one-dimensional CA with two parameters {k,r), where 
fc = |S| is the number of states, and r is the neighbourhood radius, hence EGA 
domain is defined by parameters (2, 1). There are E" different neighbourhoods 
(where n = 2r + 1) and k'' distinct evolution rules. The evolutions in this 
paper have periodic boundary conditions. 

1.2 Cellular automata with memory 

Conventional CA are ahistoric (memoryless) : i.e., the new state of a cell depends 
on the neighbourhood configuration solely at the preceding time step of ip. CA 
with memory can be considered as an extension of the standard framework of 
CA where every cell Xi is allowed to remember some period of its previous 
evolution. Basically memory is based on the state and history of the system, 
thus we design a memory function 0, as follows: 

(3) 

such that T < t determines the backwards degree of memory and each cell Si € S 
is a function of the series of states in cell Xi up to time-step t — t. Finally to 
execute the evolution we apply the original rule again as follows: 

</3(...,S-_i,S*,S-+i,...)^X*+\ 

In CA with memory, while the mapping Lp remains unaltered, a historic 
memory of past iterations is retained by featuring each cell as a summary of its 
previous states; therefore cells canalize memory to the map p. As an example, 
we can take the memory function (/> as a majority memory: 

4>Tna] Si (4) 

where in case of a tie given by Si = Sq in 0, we shall take the last value Xi. So 
(j^maj represents the classic majority function for three variables |24j . as follows: 

<^maj : {Xl A X2) V {X2 A Xs) V (X3 Axi) ^ X 
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on cells (x* 



and defines a temporal ring before calculating the 



next global configuration c. In case of a tie, it is allows to break it in favor of 
zero if Xr-i = 0, or to one whether Xr-i = 1- The representation of a EGA 
with memory (given previously in [181 117) ) is given as follows: 



4>CARri 



(5) 



where CAR represents the decimal notation of a particular EGA and m the 
kind of memory given with a specific value of t. Thus the majority memory 
(maj) working in EGA Rule 126 checking tree cells (r = 3) of history is simply 
denoted as (j)Ri26maj:3- Figure [T] depicts in detail the memory working on EGA. 
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Figure 1: Gellular automata with memory in cells. 

Note that memory is as simple as any GA local function but sometimes the 
global behaviour produced by the local rule is totally unpredictable. 

Disclaimer 

The memory mechanism considered here differs from other GA with memory 
previously reported, often referred as higher order (in-time) GAj^ These ones, 
in most cases, explicitly alter the function cj) and incorporate memory by direcly 
determining the new configuration in terms of the configurations at previous 
time-steps. Thus, in second order in time (memory of capacity two) rules, the 
transition rule operates as: x*'^^ = $(c*,c*~^). Double memory (in transition 
rule and in cells) can be implemented as: x*+ 1 = $({s*}, {s*"i}). Particularly 
interesting is the reversible formulation: a;*^^ = ip Q x*~^; reversible GA with 
memory are studied in [5]. 

2See, for example, [35] p. 118; [13] p. 43; or class MEMO in [T] p. 7. 



4 



Some authors define rules with memory as those with dependence in ip on 
the state of the cell to be updated [331 HI] • So one-dimensional rules with no 
memory adopt the form: a;*+^ = tp{xl_i,xl^i). Memory is not here indentified 
with delay, i.e., refering cells exclusively to their state values a number of time- 
steps in the past, [55]. So, for example, the cell to be updated may be referenced 
not at t but at t — 1: a;*"*"^ = (p{xl_i,xl,xl_^_i) [TB]. Again, the mapping function 
is not extended, for example, to consider the influence of cell i at time t — I: 
a;*+i = ^{x^~^ ,x\_i,x\, x-^i) as done in [4T| . 

The use of the locution associative memory usually refers, when used in the 
CA context, to the study of configuration attractors [TTJ[Tn], which are argued 
by Wuensche [31] to constitute the network's global states contents addressable 
memory in the sense of Hopfield [T^] . 

Finally, it is not intended here to emulate human memory, i.e., the associa- 
tive, pattern matching, highly parallel function of human memory. The aim is 
just to store the past, or just a part of it, to make it work in dynamics. Thus, 
working storage might replace here the use of the term memory, avoiding the 
anthropomorphic, and rather unavoidable, connotations of the word memory 
0- 



2 The basic function: EGA Rule 126 

The local-state transition function ip corresponding to Rule 126 is represented 
as follows: 



1 if 110,101,100,011,010,001 
if 111,000 




Figure 2: Typical fractal (a) and chaotic (b) global evolution of EGA Rule 126. 
(a) initially all cells in '0' but one in state '1,' (b) evolution from random initial 
configuration with 50% of '0' and '1' states. 

Rule 126 has a chaotic global behaviour typical from Class III in Wolfram's 
classification 36 . In (p_Ri26 we can easily recognize an initial high probability 
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of alive cells, i.e. cells in state '1'; with a 75% to appear in the next time and, 
complement of only 25% to get state 0. It will be always a new alive cell iff 
^Bi2Q has one or two alive cells such that the equilibrium comes when there is 
an overpopulation condition. Figure [2] shows these cases in typical evolutions 
of Rule 126, both evolving from a single cell in state '1' (fig. [2^) and from a 
random initial configuration (fig. [2|3) where a high density of I's is evidently in 
the evolution. 

While looking on chaotic space-time configuration in fig. [2] we understand the 
diflficulty for analysing the rule's behaviour and selecting any coherent activity 
among periodic structures without special tools. 

2.1 Mean field approximation in EGA Rule 126 

This section presents a probabilistic analysis with mean field theory, in order 
to search basic properties about v?fli26 evolution space and its related chaotic 
behaviour. Such analysis will offer a better spectrum where we can start to 
explore the evolution space from more useful and specific initial conditions where 
some interesting behaviours may emerge. 

Mean field theory is a well-known technique for discovering statistical prop- 
erties of CA without analysing evolution spaces of individual rules ^23, . The 
method assumes that states in S are independent and do not correlate each 
other in the local function lpri2q- Thus we can study probabilities of states in a 
neighbourhood in terms of the probability of a single state (the state in which 
the neighbourhood evolves) , and probability of the neighbourhood is product of 
the probabilities of each cell in it. 

In this way, [22] presents an explanation of Wolfram's classes by a mixture 
of probability theory and de Bruijn diagrams, resulting a classification based on 
mean field theory curve: 

• class I: monotonic, entirely on one side of diagonal; 

• class II: horizontal tangency, never reaches diagonal; 

• class IV: horizontal plus diagonal tangency, no crossing; 

• class III: no tangencies, curve crosses diagonal. 

For the one-dimensional case, all neighbourhoods are considered as follows: 

pt+i= E ^.Wpni-ft)"-" (6) 

such that J is an index number relating each neighbourhood and X are cells 
+r- Thus n is the number of cells into every neighbourhood, 
V indicates how often state '1' occurs in X, n~v shows how often state '0' occurs 
in the neighbourhood X, pt is the probability of cell being in state '1' while qt 
is the probability of cell being in state '0' i.e., q — 1 — p. The polynomial for 
Rule 126 is defined as follows: 
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Pt+i = iptqt- (7) 

Because tfBi2(, is classified as a chaotic rule, we expect no tangencies and its 
curve must cross the identity; remembering that ^fnuQ has a 75% of probability 
to produce a state one. 




Figure 3: Mean field curve for EGA Rule 126. 

Mean field curve (fig. [s]) confirms that probability of state '1' in space-time 
configurations oi(pR.i2e is 0.75 for high densities related to big populations of I's. 
The curve demonstrates also that (pRi2& is chaotic because the curve cross the 
identity with a first fixed point at the origin / = and the non existence of any 
unstable fixed point inducing non stable regions in the evolution. Nevertheless, 
the stable fixed point is / = 0.6683, which represents a 'concentration' of 'I's 
diminishing during the automaton evolution. 

So the initial inspection indicates no evidence of complex behaviour emerging 
in iy9i^i26- Of course a deeper analysis is necessary for obtaining more features 
from a chaotic rule, so the next sections explain other techniques to study in 
particular periodic structures. 

2.2 Basins of attraction 

A basin (of attraction) field of a finite CA is the set of basins of attraction into 
which all possible states and trajectories will be organized by the local function 

The topology of a single basin of attraction may be represented by a diagram, 
the state transition graph. Thus the set of graphs composing the field specifies 
the global behaviour of the system [3^. 

Generally a basin can also recognize GA with chaotic or complex behaviour 
following previous results on attractors [32]. Thus we have that Wolfram's 
classes can be represented as a basin classification: 
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• class I: very short transients, mainly point attractors (but possibly also pe- 
riodic attractors) very high in-degree, very high leaf density (very ordered 
dynamics); 

• class II: very short transients, mainly short periodic attractors (but also 
point attractors), high in-degree, very high leaf density; 

• class IV: moderate transients, moderate-length periodic attractors, mod- 
erate in-degree, very moderate leaf density (possibly complex dynamics); 

• class III: very long transients, very long periodic attractors, low in-degree, 
low leaf density (chaotic dynamics). 




Figure 4: Periodic patterns from some basin attractors. 

The basins depicted in fig. [5] show the whole set of non- equivalent basins in 
Rule 126 from I = 2 to I = 18 {I means length of array) attractors, all they 
display not high densities from an attractor of mass one and attractors of mass 
X4j3] This way Rule 126 displays some non symmetric basins and some of them 
have long transients that induce a relation with chaotic rules. 

Particularly we can see specific cycles in fig. [4] where it is possible to find: 

(a) static configurations as still life patterns (/ = 8); 

(b) traveling configurations as gliders {I — 15); 

(c) meshes (/ = 12); 

(d) or empty universes (l = 14). 

The cycle diagrams expose only displacements to the left, and this empty 
universe evolving to the stable state is constructed all times on the first basin 
for each cycle, see fig. [5] 

This way some cycles could induce some non trivial activity in Rule 126, 
but the associated initial conditions are not generally predominant. However 

•^Basins and attractors were calculated with Discrete Dynamical System DDLab available 
from Ihttp : / /www . ddlab . org/l 
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some information is useful indeed looking periodic patterns that have a high 
frequency inside this evolution space and hence for recognizing a kind of filter 
useful to get a better view of a possible complex activity in Rule 126. 

2.3 De Bruijn diagrams 

De Bruijn diagrams [23i|30j are very adequate for describing evolution rules in 
one-dimensional CA, although originally they were used in shift-register theory 
(the treatment of sequences where their elements overlap each other). Paths in 
a de Bruijn diagram may represent chains, configurations or classes of configu- 
rations in the evolution space. 

For an one-dimensional CA of order (fc, r), the de Bruijn diagram is defined 
as a directed graph with k'^'^ vertices and fc^''+^ edges. The vertices are labeled 
with the elements of the alphabet of length 2r. An edge is directed from vertex 
i to vertex j, if and only if, the 2r — 1 final symbols of i are the same that the 
2r — 1 initial ones in j forming a neighbourhood of 2r -I- 1 states represented by 
i o j. In this case, the edge connecting i to j is labeled with if{i o j) (the value 
of the neighbourhood defined by the local function) [3T] . 

The de Bruijn diagram associated to Rule 126 is depicted in fig. |6|"^| 








Figure 6: De Bruijn diagram for the EGA Rule 126. 

Figure [6] exposes that there are two neighbourhoods evolving into and six 
neighbourhoods into 1; so the higher frequency is for state 1; indicating the 
possibility of having an injective automaton; that is, the existence of Garden of 
Eden configurations HnHSHI- Classical analysis in graph theory has been applied 
over de Bruijn diagrams for studying topics such as reversibility in other 
sense, cycles in the diagram indicate periodic constructions in the evolution 
of the automaton if the label of the cycle agrees with the sequence defined 
by its nodes, taking periodic boundary conditions. Let us take the equivalent 
construction of a de Bruijn diagram in order to describe the evolution in two 
steps of Rule 126 (having now nodes composed by sequences of four symbols); 
the cycles of this new diagram are presented in fig. [7] 

''De Bruijn and subset diagrams were calculated using NXLCAU21 designed by Harold V. 
Mcintosh; available from http: //delta. cs . cinvestav.mx/-mcintosh| 
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The extended de Bruijn diagrams [53] are useful for calculating all periodic 
sequences by the cycles defined in the diagram. These ones also show the shift 
of a sequence for a certain number of generations. Thus we can get de Bruijn 
diagrams describing periodic sequences for Rule 126. 




Figure 7: Patterns calculated from extended de Bruijn diagrams, in particular 
from cycles of order (a;, 2) (that means a x-shift in y-generations). 

Cycles inside de Bruijn diagrams can be used for obtaining regular expres- 
sions representing a periodic pattern. Figure [7] displays three patterns calculated 
as: (a) shift —3 in 2 generations representing a pattern with displacement to 
the left, (b) shift in 2 generations describing a static pattern traveling without 
displacement, and (c) shift +3 in 2 generations is exactly the symmetric pattern 
given in the first evolution. 

So we can also see in fig.[7]that it is possible to find patterns traveling in both 
directions, as gliders or mobile structures. But generally these constructions 
(strings) cannot live in combination with others structures and therefore it is 
really hard to have this kind of objects with such characteristics. Although, 
moreover Rule 126 has at least one glider! This will be explained in the next 
sections. 



2.4 Filters for recognizing dynamics in Rule 126 

Filters are a useful tool for discovering hidden order in chaotic or complex 
rules. Filters were introduced in CA studies by Wuensche who employed them 
to automatically classify cell-state transition functions, see [3D]. Also filters 
related to tiles were successfuly applied and deduced in analysing space-time 
behaviour of EGA governed by Rules 110 and 54 [lH [2DJ [H] . 

This way, we have found that Rule 126 has two types of two-dimensional 
tiles (which together work as filters over fRi26)' 



the tile ii 



1111 
1001 



and 
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• the tile t2 



-0000000- 
0111110 
1100011 
0110110 

.1111111. 



Filter ti works more significantly on configurations generated by fRi26, the 
second one is not frequently found although it is exploited when Rule 126 is 
altered with memory (as we can see in following sections). 

The application of the first filter is effective to discover gaps with little 
patterns traveling on triangles of '1' states in the evolution space. Although 
even in this case it may be unclear how a dynamics would be interpreted, a 
careful inspection on the evolution brings to light very small gliders (as still 
life), as shown in fig. |8] 




This glider emerging in Rule 126 and localized by a filter is precisely the 
periodic pattern calculated with the basin (fig. |4^) and the de Bruijn diagram 
(fig. [rjD); in particular the last one offers more information because such cycles 
allow to classify the whole phases when this glider is coded in the initial con- 
dition. The next sections demonstrate the effect of filters for recognizing an 
amazing universe evolving in this CA with memory. 

3 CA (J)Ri26mA and complex dynamics 

This section discusses both relevant aspects of classic Rule 126 (1^9^126) and Rule 
126 with memory {(l)Ri26m:r)- 

3.1 Dynamics emerging with majority memory 

As it was explained in [6l [18] a new family of evolution rules derived from classic 
EGA can be found selecting a kind of memory. 




Figure 8: Filtered space-time configuration in (pRi26- 
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Figure [9] illustrates dynamics for some values of r in 4>Ri26maj 0The space- 
time configurations are also filtered to show a raw dynamics. We found that 
large odd values of t tend to define macrocells-\ike patterns |36l[23]. Even values 
of r are responsible of a mixture of periodic and chaotic dynamics. 




Figure 9: CA with majority memory (j)Ri26maj:T where 13 values of r are evolved 
and filtered over a ring of 246 cells for 236 generations. 

On exploring systematically distinct values of t, we found that 4'Ri26maj-A 
produces an impressive and non-trivial emergence of patterns traveling and col- 
liding; so yet when the memory is working as a new function, it is possible to 
recognize some fragments inherited of the original rule ipRi2e- 

We start our simulations with a single non-quicscent cell, an example of 
this space-time configuration is provided in fig. [lO| showing the first 1,156 steps, 
where in this case the automaton needed other 12,000 steps to reach a stationary 
configuration. Filter is convenient to eliminate the non relevant information 



about gliders. In fig. 10 we can see a number of gliders, glider guns, still-life 
configurations, and a wide number of combinations of such patterns colliding and 
traveling with different velocities and densities. Consequently we can classify a 
number of periodic structures, objects and interesting reactions. 



^Evolutions of 4'Rl26mar.T were calculated with OSXLCAU21 system available in http : ^ 
[//uncomp .uwe. ac .uk/genaro/OSXCASystems .htmll 
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Figure 10: Filtered space-time configuration of (t)Ri26maj:i with 602 cells, pe- 
riodic boundaries, starting from one non-quiescent cell and running for 1,156 
steps. 
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Basic primitive gliders are displayed in fig. 11 there are still-life patterns si 
and S2, and gliders gi and 52 respectively. These structures can be ordered in 

a set g4,ai26n.aj:i = { ^ 1 , S2 , ffl , 2 } ■ 



_3 



S2 



51 



92 



Figure 11: Basic gliders in 4>Ri26maj:4- Two stationary configurations (as still 
life patterns) si and S2 respectively, and two gliders gi and g2- 
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Table 1: Properties of gliders G<},„i2emaj:4- Velocity Vg is calculated as displace- 
ment between period. The linear volume is the maximum distance between two 
cells in state '1' (diameter of the set of cells in state '1'), finally mass is the 
number of cells in state '1.' 



Figure 12 shows the three more frequent glider guns (guui, gun2, and guna) 
emerging in (/)fli26maj:4- The next three guns (fig. 12 a), (b) and (c)) are com- 



bined in basic guns synchronized by multiple reactions which preserve emission 
of gliders although some of them can get other frequencies. We can include them 
increasing the set of periodic structures S<^j,i26„,a3-4 = {si: ^2, 51,52, guuj^, gun2, 
gun3,gun4}. The last gun is presented in fig. 14 as a collision of gliders. 

Basic properties of G<j>iii2emaj-i given in table fll where also glider guns 
can be classified by frequencies of glider emission as follows: 

• small guU]^ fires a gliders gi and 52 (in turn) every five steps; 

• medium gun2 fires gi and 52 (in turn) every 26 steps; 

• large gnn^ fires five gi and 52 gliders every 100 steps. 

Frequencies are related to the number of emitted gliders (in intervals) by a 
glider gun. Hence the gun2 generates one gi and another 52 glider eight steps 



15 




gmii gun2 giiiis 




Figure 12: Samples of glider guns and their compositions in 4>Ri2Qmai-A- 



before to produce the next ones, and a guns yields five gi and five (72 gliders 
100 steps before to produce the other ones. 

3.2 Collisions between gliders 

Coding glider positions to get a desired reaction is a well-known solution for 
some related problems about complex behaviours. One of them is precisely the 



problem of self-organization (by structures fig. 13 1 |15j . In this way, we present 
how each basic glider can be produced by collisions between other different 
gliders. 




(c) 



(d) 



(f) 



Figure 13: Generating basic gliders from collisions between other ones, this 
is self-organization by structures. The following reactions are illustrated as 
follows: (a) 91+92 = si, (b) si -f ^2 = ffi, (c) 51 + = .92, (d) 82+92= gi, (e) 
92 + si= 92, and (f) 91+92 = 82- 

A little bit more complicated is to obtain glider guns by collisions. Fig- 



ure 14 a) and (b) depicts the production of a gun2 from a multiple collision of 
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^ 

^5) 




(c) 



(d) 




Figure 15: Eater reactions between gliders and still lives (a, b, and c) and 
annihilation of gliders by collisions (d and e). 

A large variety of objects can be generated by collisions between gliders in 
4'Ri26maj:A- Somc of thesc objccts are useful for designing complex dynamical 
structures. For instance, fig.[l5|a), (b) and (c) shows how to delete (as an eater 
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configuration) a single glider or a stream of gliders gi and 32- Two kinds of 
annihilations are depicted as well (see fig. 15 ^d) and (e)). 




Figure 16: Black hole patterns consuming gi and (72 gliders. 





Figure 17: Examples of soliton reactions. 



Figure 16 describes a number of black hole patterns emerging in (t>Ri26maj:4:- 
Traditionally a black hole is a Life objeclj^ absorbing any glider that comes 
close to the main body (in this case, still life patterns). Its relevance consists of 
knowing how many patterns are able of attracting gliders and consuming all of 
them forever, in these cases the gi and 52 gliders. 





(d) 



(e) 




Figure 18: Examples of other collisions. 



Figure [TT] displays soliton reactions, where colliding gliders preserve their 
original forms. As known, a soliton has a small change in their phase and 
displacement; also soliton reactions are obtained with single gliders or by a 
stream of them. 

^ You can see a large classification of Life objects in ' http: //www. conwaylif e . com/«iki/| 
[index .phpl 
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Figure 19: Table of binary, multiple and other collisions. 
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The last set of examples (fig. 18 1 presents other binary reactions and some 
multiple collisions between gi and 52 gliders. Some of them are conservative 
and others produce a new number of these structures. Finally we can take 
some collisions to exploit their dynamics and controlling glider reactions (see 
an extended relations of collisions in table 19). Thus it can be obtained a full 
number of reactions for generating desired gliders. The set of collisions is useful 
as 'raw material' in the implementation of computations on 4'Ri26maj-A- 



4 Computing in 0i?l26maj:4 

Given the large number of reactions in (j)Ri2&maj-A (see table [l9|) the rule could 
be useful for implementing collision-based computing schemes [21 [2^ . 




la A 6 



a A 5 



a A 



Figure 20: Colliding interactions deriving in logic gates. 

Figure [20| illustrates the interaction of gliders traveling, colliding one another 
and implementing a Boolean conjunction as result. Initially from previous colli- 
sions we can embed logical constructions of and and not gates from this figure, 
as follows: 



for the gate (^a A h) the implementation with 
but only with one g2 glider (see tab. 19). 



'R126maj:4: 



is in fig. 18 'b) 



for the gate (a A b) the implementation corresponds to fig. 13 'a) 



for the gate (a A ^b) the implementation is presented in fig. 18 'a); but 
only with one gi glider (see table 19). 



for one fanout gate {a ^ b = a + 
fig.fTMd). 



b) the implementation is shown in 



Indeed, here we also adopt ideas developed by Rennard in his design of 
Life computing architectures |27j . Glider gi represents value 0, two gi gliders 
together represent a value 1. Two gliders 2172 traveling in positive direction 
describe the operator and one the register. Thus the register will read false or 
TRUE if they are produced successfully. 

Figure [2l]illustrates the basic reactions required to produce a primitive com- 
putational scheme in (f>Ri2emaj-A- The following set of relations is applied (see 
table fT9|: 
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Figure 21: Example of devices working as data, operator and register respec- 
tively by gliders in ^0Hi26™aj:4- 

2gx o 2^2 = e gi ^ 2g2 = gi gi ^ 2g2 = 2g2 + gi 

2gi ^ 2g2 = gi 2gi 2g2 = 2g2 + 2gi 

SO we can represent serial reactions as: 

2gi + 2g2 — e empty word 
2gi + 2g2 = gi false 
2gi + 2g2 = 2gi true. 

and a NOT gate can be represented as: 

• FALSE+252 = TKUE+2g2, or 

• TRUE+2g2 ~ FALSE+2g2- 



4.1 Constructing formal languages since gliders collisions 

in 4>Rl2%maj-A 

To be considered as a mathematical machine, 4>Ri26maj-A should compute sets 
of formal languages [TU] ■ We consider such implementation as an easy way to 
illustrate how implement some collision-based processes in 4'Ri2(,maj-A- 

Let E be a non-empty finite alphabet. Thus a string over E is a finite 
sequence of symbols from E. A set of all strings over E of length n is de- 
noted by S". For example, if E = {0,1}, then E^ = {00,01,10,11} to E^ = 
{000,001,010,011,100,101,110,111}, and so on. This way E° = {e} for any 
alphabet S. The set of all strings over E of any length is the Kleene closure of 
E and is denoted as E*. However we can write such expression in terms of E", 
as follows: 



E* = y S'' 

neZo 



This way for a binary alphabet E* {e, 0, 1, 00, 01, 10, 11, 000, 001, . . .}. 
nally a set of strings on E is called a formal language. 



(8) 



Fi- 
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Figure 22: Constructing the formal collision-based languages in (j)m2emaj:4- 
They are E° (top), (middle), and (bottom). 

Making use of gliders in 4>Ri26maj-A wc can get any set of strings of S*. For 
example, for the formal language 1* we need to code all the initial conditions 
with reactions TRUE. For the formal language (00 + 11)* the initial condition 
will be coded as: {gigi + 2gi2gi)* . To yield an arbitrary length of strings we 
increase the number of cells in the CA. All distances between gliders must be 
preserved and the left part become periodic. 

5 Conclusions 

We have enriched EGA Rule 126 with majority memory and have demonstrated 
that by applying certain filtering procedures we can extract rich dynamics of 
gliders, guns and infer a sophisticated system of reactions between gliders. We 
have discovered how a complex dynamics emerges from a chaotic system select- 
ing an adc^quate memory. We have shown that the majority memory increases 
nominal complexity but decreases statistical complexity of patterns generated 
by the CA. By applying methods as de Brujin diagrams, cycles and graph theory. 
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we have proved that Rule 126 with memory opens a new spectrum of complex 
rules, in this case a new CA with memory: 4'Ri26marA- Finally we have demon- 
strated some capacities for computing specific logical and memory functions, 
and a future work will be constructing a universal device. 
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